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Abstract. We report the first detection of a planetary 
transit by spectroscopic measurements. We have detected 
the distortion of the stellar line profiles during a plan- 
etary transit. With the ELODIE spectrograph we took 
a sequence of high precision radial velocities of the star 
HD209458 at time of a transit of its planet. We detected 
an anomaly in the residuals of the orbit. The shape and 
the amplitude of the anomaly are modeled as a change 
of the mean stellar line profile resulting from the planet 
crossing the disk of the rotating star. The planetary or- 
bit is in the same direction as the stellar rotation. Using 
the photometric transit to constrain the timing and the 
impact parameters of the transit, we measure an angle 
a = 3.9 [°] between the orbital plane and the apparent 
equatorial plane as well as a vsini=3.75 ± 1.25 kms -1 . 
With additional constrains on the inclination of the star 
and on the statistics of the line of sight distribution, we 
can set an upper limit of 30° to the angle between the 
orbital plane and the stellar equatorial plane. 
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1. Introduction 

A Jupiter mass companion having a short period or- 
bit was recently detected for the star HD209458 by 



high-precision radial velocity surveys (Henry et al. 2000, 
Mazeh et al. 2000 ) . Luckily, the orbital plane of the planet 
is close enough to the line of sight for transits to 
occur and be detected b y photometric m easurements 
(ICharbonneau et al. 200C| , |Henry et al. 2000| ). The mea- 
surement of the photometric transit across HD209458 is 
the first independent confirmation of the reality of the 
giant planets in s hort period orbits d etected by radial ve- 
locity surveys (see Marcy et al. 1999 for a review) . It leads 
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to the first estimate of the radius of a " hot Jupiter" planet 
(or 51-Peg type planet) and it strongly constrains the or- 
bital inclination i or h- 

The crossing of a companion in front of a rotating star 
produces a change in the line profile of the stellar spec- 
trum. During its transit across the star, the companion oc- 
cults a small area of the stellar disk. If the star is rotating, 
the stellar line profiles will be distorted according to the 
location of the planet in front of the stellar disk. This phe- 
nomenon was already well known from past observations 
of eclipsing binaries. It was first detected on /3Lyrae and 



Algol systems by |Rossiter (1924)| and |McLaughlin (1924) 



Recently, it has been suggested by Schneider (2000) that a 
transit by a planet could also be detected in the line profile 
of high signal-to-noise ratio stellar spectra for stars with 
high wsini. But most of the stars with a close-in planet 
have low v sin i values and in such a case the line profile 
distortion by a transit would be less than 1% of the width 
of the lines. This small effect is extremely challenging to 
detect for individual lines, but if a multi-line approach - 
like in the cross-correlation technique - is used, the mean 
effect could be large enough to be measured. Actually, any 
slight distortion of the stellar line profile changes the radial 
velocity measured from the Doppler effect, similar to the 



effect of stellar spots on rotating stars (see Queloz 1999 



for references and details). Current high- precision radial 
velocity measurements probably offer the easiest way to 
detect a spectroscopic transit of a giant planet. 

The timing and the amplitude of the drop of the stel- 
lar luminosity observed during a planetary transit mea- 
sure the radius of the planet and its orbital inclination to 
the line of sight. The detection of a spectroscopic transit 
by radial velocity measurements provides a unique means 
to estimate the relative inclination (a) between the stel- 
lar equatorial plane projected on the line of sight (called 
hereafter the apparent equatorial plane) and the orbital 
plane, as well as the ascending node of the orbit (fl p ) on 
the apparent equatorial plane. 



From a weak friction model (Hut 1981) we find that 



the tidal effect of a short period Jupiter-mass planet on 
the star is not strong enough to force coplanarity. Com- 
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parison between the coplanarization time and the stel- 
lar circularization time indicates that the alignment time 
is 100 times longer than the circularization time. The 
stellar circularization time is of the order of a billion 



years (Elasio ct al. 1996). Usually one makes the assump- 
tion that the orbital plane is coplanar with the stellar 
equatorial plane for close-in planets. Combined with the 
usini measurement of the star, this ad-hoc assumption 
is us ed to set an upper lim it to the mass of the planet 
(e.g. Mayor fe Queloz 1995 ). The shape of the radial ve- 
locity anomaly during the transit provides a tool to test 
this hypothesis. Moreover, the coplanarity measurement 
is also a way to test the formation scenario of 51-Peg 
type planets. If the close-in planets are the outcome of 
extensive orbital migration, we may expect the orbital 
plane to be identical to the stellar equatorial plane. If 
other mechanisms such as gravitational scattering played 
a role, the coplanarity is not expected. A review of for- 
mation mechanisms of close-in planets may be found in 



Weidenschilling fc Mazzari (1996) and Lin ct al. (1999) 



The amplitude of the radial velocity anomaly stem- 
ming from the transit is strongly dependent on the star's 
wsini for a given planet radius. A transit across a star 
with high v sin i produces a larger radial velocity signa- 
ture than across a slow rotator. However it is more dif- 
ficult to measure accurate radial velocities for stars with 
high vs'mi. It requires higher signal-to-noise spectra be- 
cause the line contrast is weaker. A star like HD209458 
with v ami about 4kms _1 is a good candidate for such a 
detection. With the large wavelength domain of ELODIE 
(3000A) approximately 2000 lines are available for the 
cross-correlation thus only moderate signal-to-noise ratio 
spectra (50-100) are required. 

If we use the planet's radius derived from the pho- 
tometric transit, the vsini of the star can be estimated 
from the measurement of the spectroscopic transit. Unlike 
spectral analysis, the measurement of the v sin i provided 
by the spectroscopic transit is almost independent of the 
accurate knowledge of the amplitude of the spectral broad- 
ening mechanism intrinsic to the star. A complete descrip- 



tion of transit measurements is given in Kopal (1959) for 
eclipsing binaries and Eggenberger et al. (in prep.) for 
planetary transit cases. 

2. The measurement of the spectroscopic transit 

During the transit, on November 25th 1999, we got a 
continuous sequence of 15 high precision radial veloc- 
ity measurements with the spectrograph ELODIE on the 
193cm telescope of the Observatoire de Haute Provence 
( [Barannc ct al. 199G| ) using the simultaneous thorium 
setup. The following night we repeated the same sequence, 
but off-transit this time, in order to check for any instru- 
mental systematics possibly stemming from the relative 
low position on the horizon. For both nights the sequence 
was stopped when a value of two airmasses was reached. 




508.25 508.3 508.35 508.4 

JDB-2451000 (days) 

Fig. 1. Two sequences of radial velocity residuals for the 
star HD209458 taken at the same time during the night 
but one day apart. The data are corrected for the or- 
bital motion of the planet with the Mazeh ct al. 2000 



ephemcris. Top out of the planetary transit, the resid- 
uals agree with random error. Bottom during the planet 
transit an anomaly is detected. The duration of the photo- 
metric transit is indicated by the thick dashed line. Notice 
the good timing agreement between the beginning of the 
photometric transit and the beginning of the radial veloc- 
ity anomaly. Our best model of the radial velocity anomaly 
(see below) is superimposed on the data (solid line). The 
lcr confidence level of the anomaly model is illustrated by 
the dotted area 



The ADC (atmospheric dispersion corrector) does not cor- 
rect efficiently at higher airmass. 

As usual for ELODIE measurements, the data reduc- 
tion was made on-line at the telescope. The radial ve- 
locities have been measured by a cross-correlation tech- 
nique with our standard binary mask and Gaussian fits 
of the cross-correlation functions (or mean profiles) (see 



Barannc et al. 1996 for details). 

The residuals from the spectroscopic orbit of 
HD 209458 flMazch ct al. 2000[ ) are displayed for two se- 
lected time spans in Fig. 1. During the transit an anomaly 
is observed in the residuals. The probability to be a 
statistical effect of a random noise distribution is 10~ 4 
(x2 = 53.4). The second night with the same timing se- 
quence no significant deviation from random residuals is 
observed. Note that the usual 10 ms^ 1 long-term instru- 
mental error has not been added to the photon-noise er- 
ror since the instrumental error is negligible on this time 
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Star Rotation axis 





Table 1. Star, planet and orbit parameters of our model 



star equator 



Line of sight 




Fig. 2. Illustration of the geometry of our model. The (x — 
y) plane is called apparent equatorial plane. The y-axis 
is the line of sight. Notice the definition of the angle a 
and the ascending node O p on the apparent equatorial 
plane. The angle between the star's equatorial plane and 
the planetary plane is described by On the (x — z) plane 
the thick line represents the orbital trajectory projected 
on the stellar disk. The stellar rotation axis (dashed line) 
is seen in the plane (y — z) with its angle i. 



scale, accordingly with the 40% confidence level measured 
for the non-signal model during the off-transit night. 



3. Modeling the data 

In Fig. 2 the geometry of our model is illustrated. The 
orbital motion is set in the same direction as the stellar 
rotation. This configuration actually stems from the tran- 
sit data themselves: the radial velocity anomaly first has a 
positive bump and then a negative dip. This tells us that 
the planetary orbit and the stellar rotation share the same 
direction whatever the geometry of the crossing may be 
(direct orbit). 

In this paper we decided to restrict our analysis to the 
measurement of three parameters: the angle a between 
the orbital plane and the apparent equatorial plane, the 
v sin i of the star and the ascending node £l p . The ascend- 
ing node is taken positive in the direction of the star's rota- 
tion and equal to 90° when crossing the line of sight axis. 
The timing of the transit and the impact parameter A 
(shorter projected distance between the transit trajectory 
and the star center) are set by the results of the photo- 
metric transit. The Hipparcos measurements of the transit 



(Robichon & Arenou 200C) constraint a very accurate or- 
bital period of the system. Combined with the mid-time 



of the transit by Mazeh et al. 2000, we know the observed 



transit mid-time with an accuracy of 5 minutes. More ra- 
dial velocity data at high precision would be required to 
make a full adjustment of the transit (timing and geome- 
try) . In Table 1 we list the results from the spectroscopic 



fixed parameters 

P 

a 

e 

J c 

R> 

^planet 

A* 



3.524739 ± 1.4 • 10" J d 
0.047 ± 0.001 AU 


2451508.368 ± 0.0032 
1.2(±O.1)R 
1.40(±0.17) Rj 
0.569 ± 0.004 R* 



Our best solution 



v sin % 
a 



for a = 3.9 
for a = 22 
for a — —25 



3.75 ±1.25 kms 
±3.9t£ 


100 
81 



by 



Orbit and transit data are from Mazeh et al. (2000) excepted 
for the period which is from Robichon & Arenou (2000) 
(°) computed with P of Robichon & Arenou (2000) and T,_ 
Mazeh et al. (2000) 

(*) derived from the orbital inclination angle 86.1° ± 1.6° 
( f ) ±180 [°] undefined 



orbit and the photometric transit that were used to con- 
strain our model adjustments. 

In our model we consider a spherical star in uniform 
rotation. The star is divided into 90000 cells. A Gaussian 
shape cross-correlation model with ao = S.Okms^width 
is used to model the mean individual spectral line of each 
cell where the center of the Gaussian is equal to the mean 
radial velocity of the cell. The effect of the <jq value on 
the v&mi estimate is negligible: with <ro = 2.5 kms -1 we 
would only increase the v sin i measurement by . 1 km s _ 1 . 

The integration is made by summation of the cells free 
of the planet along the line of sight. A linear limb dark- 
ening weighting (B^ = 1 — e(l — /i)) with e = 0.6 is used 
in the sum. The planetary orbit is circular. We divide the 
transit into 50 phase steps for computing the radial veloc- 
ity anomaly of the spectroscopic orbit residuals. 

To illustrate the effect of the transit geometry on the 
radial velocity anomaly during the transit, we display in 
Fig. 3 the anomaly expected for three geometric configu- 
rations and two different vsmi values. We see that the 
amplitude of the anomaly is driven by the v sin i value 
of the star. The radial-velocity symmetry of the anomaly 
of the spectroscopic orbit residuals (from the mid-transit) 
for A > impact parameter depends on the angle a be- 
tween the stellar apparent equatorial plane and the orbital 
plane. Interestingly enough, the a parameter is not tied to 
the i>sini measurement. The two parameters are uncor- 
related. In Fig. 3 we also see that the mirror trajectories 
along the apparent equatorial plane make similar radial 
velocity anomalies because fl p is ±180 [°] undefined. Nev- 
ertheless, the value of O p can be used to make the distinc- 
tion between cases a and b from Fig. 3. 
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Fig. 3. Illustration of the effect of the geometry of the 
transit on the radial velocity anomaly for HD209458: Up- 
per panel. Geometry of the crossing for 3 cases (+3 sym- 
metrical cases): (a) a = 20 [°] and VL p e [0; 90] U [180; 270]; 
(b) a = 20 [°] and tt p G [90; 180] U [270; 360]; (c) a = 3.9 [°] 
(horizontal trajectory). Lower Panel. Radial velocity 
anomalies for the geometries displayed above for (top) 
v sini=3kms _1 and (bottom) vsini=4.75kms~ 1 



The comparison of the transit model with the data is 
made by x 2 statistics with the velocity offset, the star's 
usini, and a as free parameters. To distinguish the ge- 
ometry of the two trajectories that have similar a angle 
but different fi p value, we arbitrarily set a < when Q p 
value is between and 90 (or 180-270) and a > when 
fi p value is between 90 and 180 (or 270-360). These are 
respectively the a and b trajectories illustrated on Fig. 3. 

The best set of a-v sin i solutions is listed in Table. 1 
with la confidence levels. Following our definition of a, 
any angle with 3.9 < a < —3.9 [ ] is impossible given our 
impact parameter A = 0.569 R*. Our best fit is reached 
when a = ±3.9 and il p — 0. It corresponds to a transit 
trajectory parallel to the apparent equatorial plane. 



The uncertainties on our best solution listed in Ta- 
ble. 1 do not include systematics stemming from errors 
in fixed parameters of the model. A (±0.17 Rj) change of 
the planet radius would lead to a change of the v sin i by 
±0.75kms~ 1 . Actually if the planet has indeed a larger 
radius we overestimate the vsmi value. The uncertainty 
on the star radius has a weaker effect on the v sin i mea- 
surement than the uncertainty on the planet radius. A 
change of the star radius by ±0.1 Rq would make the 
v sin i change by . 3 km s ~~ 1 only. 

4. Discussion 

We have successfully demonstrated the detection of a plan- 
etary transit using a time sequence of stellar spectra when 
the stellar line profile is distorted by the crossing of the 
planet and changes the radial velocity measurement of the 
star. We have detected an anomaly in the residuals of the 
radial velocity orbit of HD 209458 at the time of the tran- 
sit. This anomaly has been modeled with high confidence 
as the effect of the planetary transit. The data suggest 
an orbital trajectory parallel to the apparent equatorial 
plane. 

The measurement of an a different from zero and an 
ascending node off the plane of the sky (tt p ^ 0) would 
indicate deviation from the coplanarity between the or- 
bital plane and the stellar equatorial plane. In our case 
we can only argue that no evidence for non-coplanarity is 
found. However additional arguments can set some limits 
on the coplanarity level. From the statistics only, it would 
be " bad luck" to observe a non-coplanar system in a con- 
figuration such as a rj and Q p — 0. This would only 
happen for a very small set of system orientations. We 
computed the probability of finding a smaller than a max 
for a system with an angle tp between its orbital plane 
and the equatorial plane (see on Fig. 2 for an illustration 
of the ip angle). First, we have determined the relation 
between a and i for various ip configurations. Simply said 
we have computed for each configuration all the different 
spectroscopic signals that extra-terrestrial observers look- 
ing at the system from different points of view could see 
the same photometric transit. Then we have generated se- 
ries of random numbers i by recalling that the probability 
of seeing a system with a small i is smaller than the prob- 
ability of seeing a system with a large i. For each of these 
values we have used the relation between a and i to calcu- 
late the corresponding statistical distribution of a. Finally 
the cumulative distribution of a is shown on Fig. 4. With 
la confidence level we can rule out a ip angle larger than 
30 degrees. 

To get a complete description of the coplanarity be- 
tween the stellar equatorial plane and the orbital plane, 
measurement of the angle i is required. Ideally i can be es- 
timated from the rotation rate of the s tar. HD209458 has 
a quiet chromosphere (Rhk = —4.9 by Henry et al. 20()()| ) 
but significant spectral line rotation broadening (v sini) is 
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Fig. 4. Probability to measure an angle a for various con- 
figurations. The labels of curves describes the angle ?/> (°) 
between the equatorial plane and the orbital plane 
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detected. The spectroscopic usini measurements have a 
weighted mean value of 3.7 ± 0.5 km s _1 (measurements 
by Mazeh et al. 2000 and Marcy et al. private comm.) in 
agreement with the value measured in this work. The low 
chromospheric HK value suggests a rotation period of at 
least 17days (Noyes et al., 1984). With the i>sini and the 
period estimate we find i>60°. 

The statistical result on the geometry distribution of 
the system added to the constraint on i provides two 
arguments refuting strong non-coplanarity for this sys- 
tem. Better measurements of the spectroscopic transit are 
needed to get stronger constrains. Simulations show that 
for a coplanar system with the same radial velocity accu- 
racy but 4 times more data spread over the whole transit 
duration, the error bars on the measurement of a would 
be small enough to conclude that ip < 10 [°] with la con- 
fidence level. 
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